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WIATION AND SPACE TECHNOLOGY 


UDC 531.55:531.2 
MINIMUM-FUEL OPTIMAL MULTIPULSE APSIDAL INTERORBITAL FLIGHTS 


907FO0286A Leningrad VESTNIK LENINGRADSKOGO UNIVERSITETA: MATEMATIKA, 
MEKHANIKA, ASTRONOMTYA in Russian No 1, Jan 90 (manuscript received 
6 Apr 89) pp 63-67 


fArticle by S. N. kirpichnikov] 


{Text} The problem of finding minimum-fue] optimal pulsed flight 
trajectories between Keplerian orbits in a central Newtonian 
gravitational field is studied. The time of motion along the orbits is 
free. Unlike the similar problem of interorbital flight, in which the 
orbit of the flight from the initial to the final orbit must be 
constructed, one must construct « flight trajectory, begun on the 
initial orbit. and which only intersects the final orbit, in the 
considered problem of interorbital flight, so that the equalizing pulse 
rate at the finish point is not achieved. The characteristic 


maneuvering speed is minimized with the following assumptions: 1) 
“fragmentation” of the pulses is not permitted, 2) the total number of 
pulses is n § N, where N is given (N = 1, 2, ...), and 3) the flow rates 


of gases during all pulses are identical. 


Let us first consider the case of coplanar boundary orbits. Using 
argument [1], it is easy to prove the following statement, which permits 
one to limit oneself to coplanar postulation of the problem. 


Theorem 1. The trajectory of an minimum-fuel puised flight between 
coplanar boundary orbits lies in the plane of these orbits. 


The author [2-4] studied the coplanar problem of finding minimum-fuel 
optimal single-pulse (N = 1) flight. Using the methods of these 
articles, it is easy to prove the following statements, somewhat more 
general compared to the indicated articles, about the type of steady and 
optimal orbits of pulsed flights in the coplanar problem. 


Theorem 2. Any fixed orbit in the problem of an minimum-fuel optimal 
single-pulse flight always touches the finai orbit at the finish point. 








Corollars. <A limited fined multipulse orbit in the problem of an 
minimum-fuel optimal interorbital flight always touches the final orbit 
the finish point. 


The following four theorems are related to fixed and optimal 
single-pulse (N = 1) orbits in the coplanar problem of finding 
minimum-fuel optimal trajectories of interorbital flights. 


Theorem 2. The racial pulse may not. be fixed in any boundary orbits at 
all. 


Theorem 4. The tangential pulse can be fixed only if the lines of the 
apsides of the boundary orbits coincide (specifically, if at least one 
of them is circular), and single-pulse flight is apsidal, i.e., both 
apsides of the Keplerian orbit of the flight. coincide with the apsides 
of the boundary orbits. 


Theorem 5. Let the boundary orbits have a common line of apsides or at 
least let one of them be circular. The apsidal maneuvers and only these 
maneuvers are fixed single-pulse maneuvers in the considered extreme 
problem. 


Theorem 6. Let the difference of the angular distances of the 
pericenters of the boundary orbits be a multiple of 7 or at least one of 
them is circular. If the final orbit is located entirely outside 
(inside) the initial orbit, the global optimal maneuver of a 
single-pulse flight is an apsidal maneuver, which ends at the pericenter 
(apocenter) of the final orbit. 


A system of necessary conditions, derived in [2], is used in the proof 
of these theorems and the case of intersecting boundary orbits, trivial 
in the problem of an interorbital flight with free time, is eliminated 
from the consideration here and below. The following easily proved 
theorem guarantees the absence of abnormal solutions (gradients of the 
left sides of the coupling equations in the system of necessary 
conditions are Jinearly dependent) in the coplanar problem of an 
minimum-fue! optimal pulse flight between nonintersecting boundary 


mrbits. 


Theorem 7. Abnormal solutions are possible in any extreme problem of 
constructing interorbital single-pulse flights or two-pulse flights with 
free time only upon touching the boundary orbits. These solutions are 
characterized by touching all! three orbits (the initial, intermediate 
and final at a common point--the point of the start and finish). 


Ti was established above for coplanar boundary orbits with coincident 
apsidal lines that fixed and globally optimal single-pulse flights are 
apsidal. Let us study in this regard for the boundary orbits below the 
problem of the minimum-fuel optimality of apsidal multipulse 
interorbital flights. 
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Let us continue as before to assume that the problem is coplanar. It is 
easy to see that trajectories with apsidal tangential connections of al] 
Keplerian branches and the maximum cases of these connect ions-~apsida| 
tangential drifts to infinity and returns from infinity along 
parabolas--satisfy the system of necessary conditions in the case of 
boundary orbits with common line of apsides of the corresponding 
considered extreme problem. Let us also call the trajectories apsidal 
due to the brevity of them. The following statement determines the 
sufficient conditions of global optimality of apsidal trajectories. 


Theorem 8. If the orientation of one of the boundary orbits is 
arbitrary in the coplanar problem of interorbital pulsed flight (for 
example, if at. least one of them is circular), the characteristi 
maneuvering speed assumes the least value on the apsidal trajectories 
when the lines of the apsides of boundary orbits coincide. If the final 
orbit is located outside (inside) the initial orbit, the globally 
optima] maneuver with limited fiight trajectory terminates at the 
pericenter (apocenter) of the final orbit, while the difference of the 
angular distances of the pericenters of the boundary orhits should be: 
equal to (2k + 1)7, where k is an integer, for single-pulsed flight. 


Common geometric properties of transversally [5] intersecting curves on 
a plane and also some results of [i, 6-8] and methods used there to 
study the problems of interorbital flights are used in proof of this and 
all subsequent. theorems. 


Thus, let the boundary orbits now be coplanar and have a common line of 
apsides. Let us study the apsidal trajectories with N = 2. Two new 
types are possible besides the above apsidal single-pulse trajectories. 
These are primarily bielliptic apsidal trajectories. The start is made 
at the apsida of the initial orbit. (polar radius r;ji and the finish is 
made at the apsides of the final orbit. (polar radius r»}), and both 
corresponding apsida! points lie on the same semiaxis of the common line 
of apsides of the boundary orbits. The additional (second) apsidal 
tangential pulse is made on the opposite semiaxis of this line at a 
point with the polar radius p > 0. Finally, both semiellipses of the 
bielliptical apsidal trajectory change to parabolas in the final case at 
p - #, and the additional pulse becomes “zero” at infinity. The flight 
and its trajectory in this limiting case is called biparabolic. We note 
that the trajectory of biparabolic flight with such a small increase of 
the characteristic velocity can be approximated in practice by a 
bielliptical apsida] trajectory with rather large p. 


Theorem 9. Let the boundary orbits have a common line oi apsides. Only 
single-pulse and biparabolic flights can be optimal in the class of 
flights with single- and two-pulse trajectories (N = 2). It is 
necessary that r; > 4dre for optimality of a biparabolic flight. 


Tet us formulate the auxiliary obvious statement in the form of a lemma. 








Lemma. <A pulse drift maneuver, minimal according to the characteristic 
speed, with given initia] Keplerian orbit to infinity is always single- 
pulse, and the pulse is tangential and is made at tne pericenter of the 
initial orbit. 

Let us now remove the constraint N = 2. The following theorem is easily 
proved by the method of mathematical induction with regard to theorer 9 
and the lemma. 


Theorem 10. Let the boundary orbits have a common line of apsides. 
Further, let the apsides of the initial and final orbits, lying on the 
same semiaxis OL of this line, have polar radii r; and ry», respectively. 
If max{r,;/reo} < 4 upon two possible selections of semiaxis OL, some 
single-pulse flight is globally optimal in the class of minimum-fuel 
optimai flights with pulsed apsidal trajectories, and either biparabolic 
flight with start at the pericenter of the initial orbit ora 
single-pulse flight is globaily optimai in the opposite case. 
Specifically, if the initial orbit is located inside the final orbit or 


if iy < iWon (ray is the distance of the apocenter o! the initial orbit 
and von is the distance of the pericenter of the final] orbit), only 


single-pulse flights can be globally optimal. 
Let us now take the assumption of the coplanarity of the problem. 


Theorem 11. If the biparabolic flight is globally optimal in the 
coplanar problem of a minimum-fuel optima] pulse interorbital flight 
with arbitrary orientation of one of the boundary «rbits, biparabolic 
flight in the plane of the initia! orbit will also be globally optimal 
in any problem of minimum-fuel optimal pulsed flight with arbitrarily 
given arrangement of boundary orbits in space of the same dimensions and 
shape. 


The statements of theorems 10 and 11 permit considerable amplification 
for the case of circular boundary orbits. 


Theorem 12. \ Hohmann single-pulse flight (see [2]) is globally optimal 
in the class of minimum-fuel pulsed optimal flights between circular 
both coplanar and noncoplanar boundary orbits of radii r; and ro, if 
r;/ro © 4.828 ..., and biparabolic flight is globally optimal if r;/r» > 
4.828.... Any two-pulse apsidal bielliptic orbit with p > r; is more 
optimal than a Hohmann single-pulse flight for r;/re > 5.879.... 
Finally, there is a critical value Pe » ry for anv ry, re, ri/ro € 


(4.828..., 5.879...), beginning with which a bielliptical apsidal 
two-pulse flight is more optimal than a Hohmann single-pulse flight for 
any p > p,, and the economy of the characteristic speed is greater, the 


greater », p > py. All the considered maneuvers are made in the plane 


of the initial orbit. 














Comment. 


(f2 - 1)72 - 1. 


equat 


we note 


ion 


The precise value of constant 


1.828... is equal to 


Constant 5.879... is the greatest real root of the 


x4#— 9x? + 15x? + 19x + 4=0. 


in conclusion that constants 1.828... and 5.879... under 


conditions of theorem 12 play the same role as constants 11.938... and 


15.581... 


flights (see, for example, [9], section 
the conditions of theorem 12 are considerably nonsvmmetrical in the 
considered theory with respect to radii 


final] 


ho 


o>) 


—~) 


9. 


orbits. 
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in the theory of minimum-fuel optimal pulsed interorbital 


-2), but, unlike the Jatter, 


r; and re of the initial and 
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AVIATION AND SPACE TECHNOLOGY 


UDC 531.36:521.1 


ON ONE FORM OF DIFFERENTIAL EQUATIONS OF PERTURBED MOTION OF GRAVITY- 
ORTENTED SOLID 


907FO0286B Leningrad VESTNIK LENINGRADSKOGO UNIVERSTTETA: MATEMATIKA, 
MEKHANIKA, ASTRONOMTYA in Russian No 1], Jan 90 (manuscript received 
8 Jun 89) pp 71-75 


fArt cle by A. A. Tikhonov ] 


[Text] Difficulties frequently arise in construction of the transformed 
system of differential equations of motion when solving nonlinear 
problems of mechanics using approximate analytical methods, which 
provide transition to new variables. These difficulties can sometimes 
he overcome to a significant degree by using the method of variation of 
canonical! random constants. 


This method is a convenient means of constructing differential equations 
of the perturbed motion of mechanical systems, the state of which can be 
described py canonical equations. Thus, the method of variation of 
canonical random constants in combination with averaging the 
differential equations was used successfully in [1, 2] to study 
quasi-linear vibrations of a solid exposed to gravitational forces in 
the neighborhood of the position of stable equilibrium. 


4 form of notation of differentia] equations of perturbed motion of a 
gravity-oriented solid, convenient for study and which generalizes the 
corresponding equations used in [1, 2] for the case of potential] 
perturbing forces, is found in this paper. The rotational moment of the 
solid with respect to its center of mass, moving in a circular Keplerian 
orbit. in a central gravitational field, is studied on the assumption 
that. nonpotential forces that allow quadratic approximation are acting 
on the body in a central gravitational field. 


Two coordinate systems with common origin at the center of mass of the 
body are introduced into consideration: cf#(--an orbital coordinate 
system (OSK) (axis ¢ is directed along the radius vector of the center 
of mass of the body with respect to the attracting center, axis ¢ is 
directed along the tangent to the orbit in the direction of motion, and 
axis # is directed along the normal to the orbital plane, cxyz is 
rigidly bound to the body so that its axes are directed along the main 








central axis of inertia of the body, =, 4, ¢* and i, 3, k are the unit 

vectors of axes €, y. and ¢ and %, vy, Z, respectively, which satisfy the 
. . . . a b+ | 2 Pr au» ~_ ww $ > a 

equalities §=2a,i+.j+a,k, %=3,i+8.j+%3k, &=vy,i+7.J + 13K, . where 


Q _—__ ) . . 
ai, Par Th (i= 1, 3) are direction cosines. 


The following notations are adopted in the paper: « = ws + « is the 
absolute angular velocity of the body, ws = wy. is the angular velocity 
of the OSk, «° = pi + qj + rk is the angular velocity of the body with 


respect to the OSK, and A, B, and C are the principal central moments of 
inertia of the body. 


Differential Euler equations of rotational motion of the bods due to the 
effect of gravitational moment M. and of a general perturbing 


nonpotential moment Ms(|M,|<|M;,|) have the form 


d 
a (Aw)+ x (Aw) =M, + M,, (i) 


where A=diag (A, B, C), M, = 30209 X (-\G,), M, 18 the analytical! vector 


function of variables ai, Bi, yi, P, 9,7 (t=1, 3), generally non-autonomous. 





Equation (1) together with kinematic Poisson equations 





form a closed differentia! system. 


If there are no perturbing moment and fulfillment of the conditions of 
gravitational orientation B > A > C [2], the body has a stable position 
of equilibrium a; = Jo = 73 = 1 in the OSK, taken as unperturbed motion. 
It is determined by equalities vr = 7 = @ = 0 in “aircraft” angles of 
orientation (see figure). 


Posing the problem of finding a system of differential equations of 
perturbed motion in a form convenient. for study of quasi-/j inear 
Vibrations of a gravity-oriented solid, let us consider the nonlinear 
component of the gravitational moment as perturbation along with 
perturbing moment Mp: 











Following [1, 2], let us construct a perturbed differential system using 
canonical variables. Let us introduce the following canonical 


ond , aL . 
variables: 9, 9 >), p, ps, Where pi aly — Aw, sin§-+ Bw, sin g cos 64 


oL 
~ Cw, cos 9 cos 8, pr = = Aeg, P3 => ry + Boy = Bw, cos ¢— Cw, sin 9 + 
| . 
+ Boy; L = + (Aw? +Bo? +Co?) — 5 02/(A—Cyit+(B—C)%] ** ® 
Lagrange function, in which #y, @¢y, #2) 7%), and 72 are expressed by 
+, 0, 6, b, o, 6 on the basis of kinematic formulas: p = o—v¢sin 4, 
g=6cos o + 4sin¢ cos 8, r= cos gcos §—4 sin ®, 11=— sin 9, y,==sin ¢ cos 8. 


For the projections of relative angular velocity w’, we will have 


p= A™'p, — sin > cos 6, g=B"|(p, + Bw) cos o +-( 

: 4 >? + (P; +p, sin 8) X 
x sin 9 cos~' 6] — wo (sin } sin ¢ sin 8+ cos cos), r= C-! (— (p,+Bu,) x 
X sin 9 + (P, +P; sin 4) cos » cos“! 6] — w, (sin ) cos @ sin 6 — cos ) sin ¢). 


(3) 


Turning to the new independent variable u = vst, let us introduce the 
dimensionless variables 


y= (Co) p, + (1— BC) 9, 9’ =(Avy)™! p, — 4, 8 = (Buy) p,, (4) 


Formulas (4). which coincide with linearized expressions (3), divided by 
wo, are taken as the transition formulas from generalized pulses p, to 


new variables y’, g’, 6’ of the same order of smallness. Expressing 





\ ® 





z — —— ~ ~) n/ . _ : ; : - , 
Wx, Wy, Oz, @i, Bi, yi (i= 1, 3) DY YW. @»9 and retaining terms no higher than 


second order of smallness with respect to the introduced dimensionless 
variables, let us transform system (1) and (2) to the form 


(9) =F, ob +e 19/8 + (1 — de) Wo + (21 — 1) 8, 
Be (V) + OY =F, o =o +99, ia 
(VY + 2°18 (9') — dete WH) SF, SV + (1) (He —P) +H, 


Where Fe"! [42(e —8) 9 —e(1 —3 +2) 9’ + (e —8) (8 + e—1) gb’ — 24/6") — 


—(1—B +2) 9 9’+( (8—1) Y0]+M,, F,=3" [3 (e—1) 9—0'9’ + (1 —8 +2) e9/+ 


+(1—2)(¥9 — oy’) + M,, Fyset((l—t+ee’ —(’—1)(¥+3 90) + 
+ (e—1) 69" — 348'} 4+ M,, M,=(Aws)'M,,, M,=(8A05)"' M,,, M,= 


2)-4 
= (¢Awo) M,,are expressions, calculated with regard to quadratic terms 


and 6= BA", e=CA", 

System (5) at M, = O permits a canonical form of notation with 
Hamiltonian H = “tk + H,, where A, = + (A9” + B67 + Cy”) +7 + 14(B— —C) o*+ 

+ 3(A—C)#+(B—A)y?] and H, = [2 (C--B)(—9) e+ 2(A—C) o’yo 4 

+ 2 Ay’ 90+ BO’y?] + (C — B) 0 are components of second and third orders 
of smallness, respectively. Linear approximation of this system, 


corresponding to Hamiltonian Hs, i.e., an unperturbed differential 
system, has a general solution of type 


9 
y = NV; 'y. j - 4 —_ eo | - ‘ _— —1 . 6 ) 
—— uf V x; st, 9= > 3V X, cost,, 8=V2H 'x3 Rs sin t,, 
= i=-] 


2 9 
t =- \’ ; o e or : ara * ny - 
? pa k,2,,.V x, cos noe > k,32,V x, sint, 6’ = V2 Bx, cos *,, (7) 


i=mt 


where tW=;(u + v5) (j=1, 3); x;, vj; (j=1, 3) are canonical arbitrary 


constants and k, (j = 1, 3) are unperturbed natural frequencies of the 
body: 














rx Va27VaA—d, h=V3U—e0, a=143(8—2) +0/4, 
ba 4e%(8—1)(8—e); 0,2 V2 xk VAL, oy =—V2K'V AT 


(= 1 Be (Ap A) y(t 8 +8) (i=1, 2) 
EY ed ot) ee ed es | 








The inequalities 4,>0 (i=1, 2), m>0, ‘x2 <0, o1>0, o12<0, 021<0, o22<0. 


are valid in the region B > A > C. 


Solution of system (5) by the method of variation of arbitrary constants 


xi, Vi (i = 1, 3) is sought in the form of (6), where xj and vj are new 
unknown functions subject to determination. 


Observation. It follows from the kinematic equations of system (5), 
written in variables x; and v,, that the values of x’, and v’, have the 


order O(¥<3,) and O(/x,), respectively. Taking these estimates into 
account and also taking into account that the differential equations of 
perturbed motion with respect to conditions of postulation of the 
problem should contain terms no higher than second order of smallness 


with respect to ¢, ¥, 6, 0’, o, & in derivatives (4’)’o, (9’)’6, (6) o, 


contained on the left sides of the dynamic equations of system (5), and 
one should accordirzi: leave only the following terms of second order of 


sma] lness: 9 yi K Sy V x; sin "Oo > Riz V x,cost,; — oksV 2B ' x; sin 75, 
iol i=l 


where ~ and #¢ are determined by formulas (6). 


Turning in system (5) to variables x; and vj; according to transform (6) 
and (7) and solving it with respect. to the derivatives, we find the 
desired system of differential equations of perturbed motion of the 
body, which can be written in the form 





’ 6H > -147 . —l ay 
x; = -= + 2 (— 1) Vx; (oy, 31 94 'M, SIN 7, — C2, 3-192 'M, cos AP 
’ 0H. =e . -] a ~l ar ° 
v= ae + (h, V x;) ‘(- 1)‘ [o), s19 M, Cos t, — 52, 3-32 M, sin =], 
1 OM; V2Bx.M ‘Oy VB . 
+3=— Fy, + 2Bx,M, COS Ts, 03 = “Oxs kV 2x, M, S1N 73. 


for H, in variables «;, vj has the form 


10 








H, = (A,X, + 2X2) V x, cos t + Asx, V x; cos (24 + 13) + hyx,V x; X 
X C08 (2% — ty) -+ Asx, V x, COS (2 ty %) + Ay x, V x3 COS (2 t9 — ts) + 
+ hyV X,XoX; COS (t, ++ t; — %) ++ hgV X,XoX3 COS (t, + t3 — %)) + 
+ AV x,%_X; COS (4 bt, —%5) + Ayo V X X24, Cos (=, + % + %), 
where ha=VB/2 [(e8™t — 1) (bet1a — G22) O20 + S4e/2] (a =1, 2), 
hg =V B/8 [(e8-—-1) Sn, (R; 94,—42,+(— 1)? ks'3,,) — 01;/2+(— 1)P3* RX 
Xx ks '0,, [(1—e)¢,, — k,32,)] ({=1 ana B=3, 4; 2=2 for 8 =5, 6), 


anaael 2 
VBR] oo 1) o,—(—1)ta5 "8" »(- 1)'R3—191, 3-1[(1 —£) 04) —Ry2,] + 
im} 
+ (— 1)tS! (1 — €) (35999) — 941922) + 944512] (7 = 7, 8), 
h,=V Bis [(e3* — 1) 35 + (—1)*RF 18" [(Ry Re 8 — €) (94:83 + 949921) — 
— (1 —¢) (2, + 22) 6,,9;2] — 711719] (v==9, 10), 


93 = R44 4Fy9 1 RyFy2Fq1 — 2 F949 29. 


The convenience of using the derived differential equations of perturbed 
motion in the form of (8) for study of quasi-linear vibrations of a 
solid is determined by the following circumstances: 1) the structure of 
these equations is such that their right sides are divided into two 
groups of terms, from which the first, written in canonical form, 
corresponds to “canonical” perturbations, caused in the given case of 
nonlinear quadratic terms at the gravitational moment, while the second 
group corresponds to "“noncanonica]” perturbations due to moment Mp? 


generally being of an arbitrary nature; and 2) system (8) has a form 
that permits direct use without laborious computations of asymptotic 
methods of nonlinear mechanics for study of it. 


Equations (8) are canonical in the absence of a perturbing moment ‘My and 


coincide with the widely used equations in canonical variations [2]. 
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NUCLEAR ENERGY 


UDC [621.311.25:621.039]:620.193.4.091.45 
Some Features of Corrosion Processes in Chemical Cleaning of AES Heat Exchanger Equipment 
907F0210A Moscow ELEKTRICHESKIYE STANTSII in Russian, No 2 Feb 90, pp 13-16 


[Article by N.N. Berezkin, engineer, Yuzhtekhenergo; and S.S. Kirilyuk, C.Chem.S., Lvov State Medical 
Institute] 


[Text] An aqueous solution of a synthetic low-molecular organic:acid (NMK) is often used as a cleaning 
agent during the process of chemical cleaning of carbonaceous deposits on heat exchanger equipment. Use 
of this acid in AESs involves the particular characteristics of the structural materials in the primary 
equipme nt»? As was shown in ”4, steps need to be taken to reduce the rate of corrosion of the 
equipmen: during cleaning in connection with the aggressiveness of the cleaning solution as compared to 
the metal from which the deposits are being removed. 


One method of reducing corrosion in the copper alloys used as the primary material in the manufacture of 
heat exchangers is the protective screen method. Before chemical cleaning, protectors — steel plates — are 
mounted along the loop being washed. According to the data in the literature one may say that the rate 
of corrosion of the protectors increases sharply, taking most of the load of the corrosion process upon 
themselves, and because of this the rate of corrosion in the primary material of the heat exchanger, in 
particular MNZh-5-1, should be reduced considerably; MNZh-5-1 is used mainly in the manufacture of the 
tubing systems for the condensers in turbines for high power generating units. 


However, we should note that as the steel protectors dissolve, an additional significant amount of 
trivalent iron ions land on the loop being cleaned, and the effect on the nature of the corrosion process is 
clear: there is a dramatic increase in copper alloy corrosioi, i.e. intense wear of the heat exchanger during 
chemical cleaning. Even if protective screening is not used, steel structures along the loop — condenser 
tube plates, tubing intakes and drain pipes — may take the place of the protectors. 


As experience has shown in the use of periodic chemical cleaning, a lot of cleaning solution (and, hence, 
money) is lost at the ends of condenser intakes and drain pipes and in the places where the tubes roll 
around in the tube plates, i.e. where intense local corrosional wear takes place. 


With the idea of reducing the impact of trivalent iron ions on the corrosion rate of copper alloys in acidic 
cleaning solutions, we recommend the introduction of compounds that facilitate the transfer of trivalent 
iron to the divalent state *. Thiourea is widely employed for this purpose, in addition to a number of 
other compounds. However, the corrosion process in which steel and MNZh-5-1 alloy are present together 
in acidic solutions with admixed reducing compounds inhibitory to acidic corrosion has not been studied 
enough. There is cause for interest in carrying out an investigation in this area. 
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Figure 1. A Graph of the Corrosion Processes For MNZh-5-1 and Steel 3 Being Present Together in a 7% 
NMK Solution, 7’ = 333 K: 1— MNZh-5-1 corrosion rate (No Steel 3 present); 2— MNZh-5-1 corrosio. 
rate; 3— effect of thiourea (0.25 g/dm®) on the MNZh-5-1 corrosion rate; 4— 6 — protective action of 
acidic corrosion inhibitors (US-1, R-12, R-14, R-15) on MNZh-5-1; 7— Steel 3 corrosion rate (no MNZh-5- 
| present); 8 — Steel 3 corrosion rate; 9— effect of thiourea (0.25 g/dm*) on the Steel 3 corrosion rate; 
10 — 13 — protective action of acidic corrosion inhibitors (US-1, R-12, R-14, R-15) on Steel 3 


Figure 2. A Graph of the Corrosion Processes for MNZh-5-1 and Steel 3 Being Present Together in a 7% 
NMK Solution, 7' = 313 K: Notation the same as for Figure 1. 


The methodology for conducting a laboratory experiment has already been described 152° Steel 3 and 
MNZh-5-1 samples are simultaneously placed in chemical beakers containing a 7% NMK solution, to 
which thiourea is added (0.25, 0.5 and 1.0 g/dm?, tests 3 — 5, see Table), as well as certain acidic 
corrosion inhibitors that we synthesized ourselves (US-1, R-12, R-14, R-15), having a constant 
concentration of 0.5 g/dm?® (tests 6 — 9, see Table). In addition, Steel 3 and MNZh-1-5 were placed in 
beakers with NMK separately (test 1) and together (test 2) without adding thiourea or inhibitor in order 
to study the effect of MNZh-5-1 on the corrosion rate of Steel 3 and vice versa. 


The Table shows the results of the corrosion tests at 298 K, while Figs. 1 and 2 show them for 313 and 
333 K. 


As follows from the obtained data, the corrosion rate of Steel 3 in a 7% NMK solution with no additives 
or MNZh-5-1 is practically invariant over time. By contrast there is even a small decrease in the 
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oTeeL 2 rath corrosion rate of the MNZh-5-1 alloy in the same 
| on situation, which agrees with the earlier work °. 























rest 7 of ” 2. % ” of | 2% rene h The test were conducted under static conditions 
19/m*h mh without stirring of the acid solution. The Steel 3 
0.09 0 0.09 r 24 and MNZh-5-1 corrosion rates were found by 
0,08 0 0,08 0) 48 averaging five parallel tests, as is instructed in 
aa : Y : % the CEMA 4815-84 standard. 
0,09 0 0.07 0 120 
2° pis ; ye 0 24 When Steel 3 and MNZh-5-1 are both present, 
016 0 yo : i the corrosion process changes. First of all, the 
0,16 0 0,09 0 120 MNZh-5-1 corrosion rate “stabilizes” at a higher 
3 ye . Pit a a level, and secondly, there is observed a 
0.11 99 0.12 —64 72 substantial increase in the Steel 3 corrosion rate 
0,11 29 0,14 —60 96 (by a factor of 1.6 to 2). In addition, there is a 
4 ri he 4 yo a * + significant amount of reduced copper noted on 
0.18 — 0.09 7 4s the steel plates, which is difficult to remove by 
mechanical means. 
0,11 32 0,03 63 0 
5 . , ye s Fy Introducing thiourea into the NMK solution 
001 95 0:02 83 .. turns out to have a differing effect on the 
te ro ‘ie *. b. corrosion process depending on its “dnt the 
' ' For a thiourea concentration of | g/dm” the 
. os 5 yt ro 1 degree of protection for Steel 3 (the decrease in 
ye 79 0,01 86 72 the corrosion rate) fluctuates between 47 and 
yt a yt os b 8.5% by the end of the experiment, while the 
7 0,02 83 0,005 94 24 degree of protection is chiefly negative relative to 
yy a yee : . MNZh-5-1, i.e. there is an increase by a factor of 
0,03 78 0,009 90 96 1.6 to 1.9 in the corrosion rate. The most 
‘ ry 2 yy! s i” beneficial thiourea concentration in this instance 
0.02 86 pn pnd 48 should be recognized at 0.25 g/dm? (test 10, 
0,03 8) 0,008 90 72 Table), for in this instance the Steel 3 corrosion 
yt . yr z i rate was considerably reduced. Nevertheless we 
9 0,02 84 0,0! 92 24 need to point out that in all the tests using 
ye rs 001 7 . thiourea (tests 3 —5, Table) there was significant 
0.03 78 0.01 90 96 precipitation of elementary sulfur and, in 
0,03 80 0,01 86 120 addition, to some degree or another all the steel 
10 She 4 ay vl a plates had reduced copper deposits on them. 
0,06 64 0,009 90 79 The effectiveness of chemical cleaning is reduced 
Boe 3 yt D when precipitate falls into the recessed areas of 


the loop, and the presence of cupric deposits on 
the steel surfaces leads to an increase in the 

*Steel 3 and MNZh-5-1 are present together corrosion rate of the loop as a whole and to the 
need for mechanical cleaning or other measures. 





The best idea, as we have proposed, is to apply a reliable protective inhibitor to the entire loop to be 
cleaned?. It was earlier shown”?® that a number of the inhibitors we synthesized effectively protect 
against oxidic corrosion in ferrous metals as well as in the MNZh-5-1 alloy. For reduced corrosion rates in 
the most highly stressed sections of the loop in chemical cleaning of the heat exchanger, it is relevant to 
study the inhibitory process when both Steel 3 and MNZh-5-1 are present in the NMK solution. 
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The addition of 0.5 g/dm? of some inhibitor gets us from a 78 to <1 86% degree of protection for Steel 3 
(Table, tests 6-9) and from 83 to 94% for MNZh-5-1 alloy. The inhibiting effect for higher 
concentrations of these compounds is to be found in '. We should note that there is a complete absence 
of reduced copper on the steel samples, consequently the electrochemical process of reducing copper has 
been completely suppressed. 


Negligible concentrations of thiourea lower somewhat the corrosion rate of Steel 3 in the 7% NMK 
solution. 


In order to study possible synergy, a separate series of NMK solutions were prepared that had thiourea 
and oxidic corrosion inhibitors added at the same time. As was demonstrated in the investigation (Table, 
test 10), in spite of the absence of reduced copper on the steel plates, the degree of protection of the 
control samples did not increase, i.e. no synergy was discovered (no increase in the inhibitory effect). 


In order to study the effect of cleaning solution temperature on the corrosion process, we also performed 
tests at the higher temperatures of 313 K and 333 K, which are sometimes employed when acid cleaning 
carbonaceous deposits. As follows from the data (Fig. 1, 2), the nature of the corrosion process changes 
almost not at all with increasing cleaning solution temperature in this case, i.e. the presence of thiourea 
facilitates increased corrosion wear in MNZh-5-1 alloy. 


Thus we have studied the effect of both Steel 3 and MNZh-5-1 being present in an NMK solution on the 
corrosion process. We have remarked on the role of protective screens in chemical cleaning of 
carbonaceous deposits from heat exchangers. The ineffectiveness of using thiourea under actual 
circumstances has been recognized. The effectiveness of using a reliable protective inhibitor has been 
confirmed. 
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NUCLEAR ENERGY 


UDC [621.311.25:621.039].001.572 
KIPR Time-Dependent Software Package for Investigating the Dynarnics of VVEh-1000-Equipped AES 
907F0210B Moscow ELEKTRICHESKIYE STANTSII in Russian, No 2, Feb 90, pp 16-20 


[Article by I.M. Arshavskiy, C.Tech.S, V.M Dmitriyev, D.Tech.S., S.B. Korolev, engineer, and A.Ye. 
Kroshilin, D.Tech.S.; VNIIAES] 


[Text] A great many studies have been devoted to investigating the thermal and hydraulic processes 
going on in non-stationary and emergency modes of operation of VVER-1000-equipped AES. Most of the 
mathematical models and programs that have been developed are directed toward reaffirming the 
operating safety of the power generating units 


The problems that arise during power generator operation place particular demands on a mathematical 
model. We will list a few of these. 


1. We need to have a mathematical model of the power generating unit as a whole, including the 
primary process systems, automatic reguiators and measuring devices. 


2. The mathematical model should be able to be readapted for each different actual AES power 
generating unit (taking into account the individual technical implementations, adjustments to the 
automatic regulating systems, etc.), as well be able to process output data for monitoring by AES 
technologists. 


3. The mathematical model should be both fast and accurate. 


Satisfying the latter requirement is accomplished by developing a set of mathematical models, each of 
which has its limited area of application while ensuring efficiency and accuracy. The combination of 
mathematical models allows us to have a wide range of applicability of the mathematical model as a 
whole while maintaining its speed. 


The complex of computer programs that we have developed were designed on the modular principle. All 
the AES process equipment and physical processes were nominally divided into groups, each of which was 
described by a separate functional n.vdule (subprogram). Most of the functional models had variant 
implementations differing in their detail and accuracy of describing the corresponding object or process 
and, as a result, the computational characteristics of the program: its speed and the amount of computer 
memory needed. The software package included the following functional modules: 


N—a subprogram for neutron physics calculations; H — a subprogram for modeling the thermal and 
hydrodynamic properties of the coolant in the primary loop of the AES; TEPL — a model describing the 
temperature fields in the fuel elements, reactor shell and piping; TPS — joining the models of the primary 
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Figure 1. Structural Diagram of the KIPR Software Complex. Key: 1 — Control; 2 — Data; 3 — Display; 
4 — Magnetic Disk; 5— Common Area in Computer Main Memory 


loop process system and reactoi core emergency cooling system, including the pumps for emergency 
reactor flooding and aftercooling; PG —a model of the steam generators and steam lines; 72— a model of 
the secondary loop thermal equipment, including the turbine, condenser, deaerator, steam superheater, 
etc.; and AUTO-—a functional subprogram for joining all automatic regulating systems of the primary 
and secondary loops of the power generating unit. The operation of the automatic systems is described 
using the real control laws implemented by actual regulators, their adjustment parameters, dead zones, 
servomotor inertia, etc. AWN is the model of the AES power system. 


To ensure operation of all the program modules in the context of the overall mathematical model, w: 
standardized the content and form of data exchange between program modules, and the rules governing 
their operation. All data exchanged between the functional modules is located in a common area of the 
computer’s main memory. A dedicated subprogram monitor controls the operation of the functional 
modules, which organizes calculation of the initial steady-state conditions of the power generating unit 
and calls the functional modules in cyclic ordet (modeling a process that is non-steady with respect to 
time). The monitor also provides interactive operation with the model, so that the researcher has the 
opportunity to output the model’s results in graphic and text form to any number of displays and vary 
the operating mode of ihe program in an interactive manner. In addition to the functional modules and 
control subprogram as a part of the software package, there are also a number of service programs, 
including /NJT —a module that accompanies the initial data needed by the software package. Initial data 
includes several thousand numbers, listed in tables, laden with comments and special indicators. The 
INIT service program substantially simplifies the preparation of initial data for calculations and allows us 
to construct a mathematical model for an actual power generating unit. 


The structure of control and data transfer of the software package is shown in Fig. 1. 


This organization of the software complex has the following important advantages permitting efficient use 
of this tool for operational tasks: 

eit enables us to broaden the content and, consequently, to increase the detail and accuracy of the 
model; ' 

eit gives us flexibility in the modeling system (more accurate, slower models used only in certain 
circumstances, or, for the majority of calculations, faster and simpler models); 

ethere are the service programs that allow AES technologists to readjust the mathematical model for 
actual power generating units and process the output data (including interactively). 


We will now give a more detailed description of the fundamental functional modules of the software 
package. The calculational complex consists of 4 interchangeable versions of a model of the thermal and 
hydrodynamic state of the coolant in the primary loop of a power generating unit (H/— H4). 
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1. The primary job of functional module H/ is to model the processes associated with large-scale flows in 
the primary loop of the power generating unit (including the Maximum Design Accident). What is being 
determined in this case is the correct description of the effect of pressure wave propagation. 


One-dimensional time-dependent control of the conservation of mass, momentum and energy in a 
homogeneous steam-water flow is integrated in. It uses the calculational method of ', which is a 
generalization of the Godunov numerical scheme? for the case of more complex (as compared to an ideal 
gas, for which the Godunov method is traditionally employed) control of the state of the steam-water 
mixture. 


We note the advantages of the numerical scheme H/ is based on as compared to other schemes for 
hydrodynamic calculations (e.g., the Lax-Vendroff scheme, or the method of characteristics): the H1 
numerical scheme describes more accurately the wave flow of coolant for a minimal number of 
differencing steps; conservation of mass and energy is built in to the numerical scheme, which is especially 
important for investigating loss-of-seal in a loop; and the numerical scheme is highly effective and stable 
for supersonic flows as well as those approaching sound speed. 


2. The primary job of model H2 is also a description of the processes associated with loss-of-seal in the 
primary loop. A one-velocity two-temperature model of coolant flow is employed. It is assumed that the 
ternperature of the steam is identical to the saturation temperature and it considers only the most 
significant effect — the possibility of water making the transition to a superheated state in the initial 
stages of boiling. 


The kinetic equations for steam generation include a function for the steam generation intensity, which 
follows from a generalization of the large amounts of experimental data. It considers boiling at the walls 
of the tubing as well as within the volume of the superheated liquid on solid impurities. 


Article 3 presents the calculational model on which functional module H2 is based in more detail. 


3. When modeling time-dependent AES operating modes with characteristic times on the order of tens of 
seconds, use of mathematical models H/ and H2 is involves significant outlays of machine time. This is 
in connection with the Curant limit At < Ar/|v+c], which determines when to “compress” the 
propagation of acoustic disturbances. 


We may point out that there is a rather broad class of slow time-dependent modes in which wave effects 
are not decisive and use of mathematical models H/ and H2 is not justified. 


Mathematical model 3 is based on a quasi-stationary calculation method for the coolant flow in the 
loops of a power plant 4. This method describes the quasi-stationary flow class in which variations in 
velocity and pressure over time are governed primarily by the time-dependent boundary conditions and 
variations in the external parameters. Nondimensional analysis shows that an integrated, simplified 
system of hydrodynamic equations is accurate enough for this class of flows: time-independent equations 
for velocity v and pressure p, and time-dependent equations for the density p. The Curant stability 
condition for this case has the form At < Az/|v|, which enables us to construct solutions with larger 
time steps. 


For description of the effects of boiling and condensation caused by pressure changes over time, we use a 
constant-volume idealization, appropriate to which we introduce the unknown parameter dp*/dt, which 
characterizes the speed of pressure changes over the loop. 
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For integrating the continuity equations we use a conservative numerical scheme 5 based on analytical 
solution of the decay of an arbitrary density discontinuity in a two-phase two-velocity steam-water flow. 
The fundamental steam-water flow regimes are considered: bubble, slug, eddy, disperse, and transitional 
regimes. 


Thus functional module //3 ensures conservation of mass and energy, takes into account the expansion of 
coolant because of changes in pressure over time, and correctly describes the effects of velocity 
nonuniformity with a small number of calculation steps. 


4. In the modeling, a fairly broad class of nonstationary AES modes do not require a high degree of 
accuracy in descriptions of the coolant flow in the primary loop of the power generating unit. The slow, 
nonstationary processes with characteristic times from hundreds to thousands of seconds in which the 
coolant remains one-phase are an example of this. 


We have developed a simplified hydrodynamics mathematical model H4 for description of these regimes, 
and for performing multivariate quantitative calculations. 


The primary loop of a power generating unit is represented by three elements (the mixing chambers and 
the pressurizer) with averaged parameters. Changes in the flowrate of coolant over the length of a loop 
and core channels are ignored, and we also do not take into account the time it takes for coolant to pass 
through the loops and core channels. The key to the method is the following assertion: 


From the equations of conservation of mass and energy, we obtain a relation for the rate of change in the 
pressure in the NKS and VKS [not further identified] with respect to time. 


The total system of equations for determining the pressure and flowrate in the loop includes equations of 
the form Ap = €G? (where € is a known coefficient) for the loop and core channels, the equation of 
conservation of energy for the pressurizer, and balancing equations of conservation of energy for the loop 
as a whole. We propose an implicit calculation scheme for solving these equations, which allows us to 
increase the time integration step up to 1 s. This makes the functional model quite fast. In additional, 
mathematical model H4 provides enough accuracy for engineering calculations, which follows from a 
comparative analysis of its operation, experimental data, and calculations from the models H3 and H1. 


The software package includes two interchangeable models for neutron physics processes in the reactor 
core (Nos 1 and 2). 


1. This model is designed to describe the processes in which spatial neutron physics effects are 
insignificant. Mathematical model No 1 is based on integrating the equations of point kinetics in a quasi- 
stationary approximation with delayed neutrons being broken down into six different groups. It also 
considers fast feedback on fuel temperature, water temperature and boron concentration. For a 
calculation of the reactivity induced by the control rods we use an experimental dependence of the 
integral effectiveness of a control rod bundle as a function of the depth of penetration into the core. 


2. Module No 2 is designed for a description of the spatial distribution of the neutron flux in the core. 
The core is divided into 10 calculation regions, and in each of these the neutron flux is considered to be 
piecewise-continuous. A quasi-stationary one-velocity diffusion approximation is used, with delayed 
neutrons being broken down into six different groups. 
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medium is represented in the form of a piecewise- 
continuous function that is independent of the 
spatial coordinates at the boundaries of each of the 
calculational regions). 


The very real advantage of this model is its 
efficiency, which makes it possible to use neutron 


kinetics model No 2 for multivariate calculations. 


A mathematical model of the horizontal steam 


generator has been developed in two interchangeable 


variants (functional modules PG/ and PG2), 
differing in the detail in their description of thermal 
and hydrodynamic processes, which are realized in 
steam generators in the time-dependent operating 


Figure 3. Graph of Variation in Parameters for 


Switching Off One of the Main Circulating Pumps at 


YuUAES: 


1 — Flowrate of steam from the “emergency” 


steam generator; 2 — Flowrate of feedwater into it; 3 — Level 


in the steam generator with respect to a baseline of 1000 
mm; 4 — difference in temperature between the hot side and 
cold side; 5 — pressure in the steam generators; 6 — pressure 
in the hot side: — experiment, — — — 

— — calculation for the “emergency” loop, 

—xX—x—x— — calculation for the “non-emergency” loop, 


—_ came ams as — level readings without consideration of inerti 


Figure 2. Graph of Pressure in the Second Loop (1) 

and in the First Loop (2), and Variations in 

Temperature at the Reactor Inlet (3) and Outlet (4) 

as a Function of ‘Time in the Closed SRK Turbine 

Mode (With Triggering of the Emergency Shutoff): 
— DINAMIKA program calculations; 


KIPR program calculations 


For integration of the neutron kinetics equations we 
use a method based on expansion of the neutron flux 
density function in terms of eignenfunctions of the 
Laplace operator, which allows us to reduce the 
problem of integrating the initial differential 
equations to solution of a system of linear algebraic 
equations in terms of the average neutron flux over 
the core (here the breeding characteristic curve of the 
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Figure 3. Graph of Variation in Parameters for a Np, ‘Vo 
False Triggering of the Scram System at the KaAES: 100 





1 — Neutron power; 2 and 3 — Pressure in KD and GPK; 4 





























85 
and 5 — temperature at the reactor inlet and outlet, 
— experimental, — — — — — calculation 70 
regimes of a power generating unit. Both models 551 
(PG1 and PG2) include a description of the steam 
generators and steam lines connecting them to the 40) 
turbines. We introduce a description of the more 
complicated model PG2. 25 
Prax» 
In modeling of steam generators, the most complex 6,7 


problem is determining the water level. The PG2 
distributed model was developed for a detailed 
description of the processes in the two-phase volume 61 
of the steam generator. It takes into account the 

following primary effects that substantially influence 5,6 -— 
the thermal and hydraulic state of the steam aN 
generator and the water level readings in time- 5,5 . 
dependent processes: phase slipping in the steam- 

water volume of the steam generator, natural , 
circulation of coolant, variations in pressure drops 49 
over the submerged tube plate, the complex "0 
geometry of that plate, and the presence of flow at 

the edges. 






































The presence of the submerged tube plate and natural circulation of coolant in the steam generator results 
in a number of instances in the level gauge readings differing substantially from the actual amount of 
water in the generator. In determining the level gauge readings, model PG2 calculates the drop in 
pressure which will be measured by the gauge. The equalizing vessel with pulse lines is also modeled, 
which enables us to describe the effect of needle inertia in the level readings. This ensures a fairly high 
degree of accuracy in predicting level readings during time-dependent processes. 


The mathematical model for the unit’s automatic regulation system is implemented in only one variant. 
The functional module AUTO models the work of the primary AES power generating unit regulation 
system in normal, transition and emergency modes, corresponding to complete or partial power 
reductions. The model describes the control system of the primary loop of a reactor, the pressurizer, the 
secondary loop and turbine. Analysis of the values of controlled and observed variables, realization 
appropriate to the structure of the regulation and control algorithms, and the formation and transfer of 
vector control actions to the other subprograms are all functions of the AUTO module. The AUTO 
module compares current parameters with set values and generates safety signals or interlock signals, as 
well as diagnostic information. 


The modeling of the regulation system was carried out in volume, allowing us to reproduce the operation 
of the power generating unit over the range of powers and to realize all normal and emergency modes. 
Operation of the regulators were described by differential equations appropriate to the regulation 
guidelines that shape them. Static as well as dynamic regulator parameters were taken into account. 


One of the most important stages in the development of each functional module included in the 
mathematical model software package is comparing the results of its operation with experimental data 
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and with calculations from more accurate models. This kind of independent testing lets us verify that the 
approximations made in the mathematical models are fair, and to refine the cross-coefficients in the 
correlations used. 


In order to check on the package as an adequate mathematical model of the object, we conducted a 
comparison of the calculations and experimental data obtained during dynamic tests in AES. A 
comparison was also made with the DINAMIKA program’. 


Fig. 2 shows a graph of the changes in unit parameters when the turbine SRK are closed (with triggering 
of the emergency shutoffs). The graph in Fig. 3 illustrates thy calculational and experimental results in 
which one of the four main circulation pumps that operate under nominal power is switched off, which 
happened at the YuUAES®. The graph in Fig. 4 characterizes the results of the false plunging of the 
control rod bundle into the core, as conducted at the KaAES No 2 unit. 


The good agreement between calculation results and dynamic test data (about 30 different modes were 
processed in all) opens up the possibility of using the KIPR software package for meeting a wide range of 
needs arising in the process of operating an AES, including the development of instructional methods for 
the operating personnel. 
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[Text] VK-100-90 turbines operate at p = 90 kgf/cm”, T = 500°C and at present have logged about 
300,000 h. 


In accordance with the technical reequipment plan of USSR Minenergo, all condenser turbines of this 
brand are to be written off, and turbines converted over to central-heating uses wili be modernized, in 
which the replacement of the high-pressure cylinder in the assembly is stipulated. 


This article presents material about inspection of the metal in the turbine as a whole, in both the high- 
pressure and low-pressure cylinders. 


High-Pressure Cylinder Housing. By inspecting the condition of the metal in this turbine after 200,000 h 
of operation, we have shown that the high-pressure cylinder may be operated reliably for up to 500,000 
h'»?_ This conclusion was drawn on the basis of defectoscope monitoring, studies of microstructure zones 
over the operating temperature regimes of the turbine, analysis of the fine structure, and study of the 
phase analysis and mechanical properties. 


Recently we conducted additional studies with estimates of the crack resistance of 20KhML steel and 
investigation of the condition of the turbine metal after 275,000 h. As is known, the cyclic ultimate 
strength and crack resistance criteria should be numbered along with the static and long-term strength as 
limiting criteria of the condition of the turbine housing. 


The primary method for estimating the crack resistance in the process of operating the turbine is 
magnetic powder defectoscopy (MPD) at the radius transitions. After 300,000 h of operation the housing 
parts were subjected to monitoring three and four times. Thirteen turbines were subjected to this 
monitoring, installed in three different electiic power stations. The turbines operated in a base-load 
regime. The operating period of all the turbines was practically identical. The number of starts ranged 
from 200 to 500. 


Table 1 presents the data on the number of detected cracks, and Table 2 gives the depth characteristics of 
the cracks. 


All the cracks occurred only on exterior surfaces. 
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Tabl« 1 





Number of cracks over the operating period, in thousand of hours 
























































Control object 
up to 100 100 — 150 150 — 200 200 — 250 

HPC Housing 5 6 41 35 

HPC Cover 21 20 20 18 

Stop Valve 6 2 13 ll 

Regulator Valve 19 18 25 13 

Table 2 
Crack Number of cracks over the operating 
period, in thousands of hours 
—— “— to 100 — 150 — 200 — 
inmm | “foo 150 200 250 

HPC Housing up to 5 ~ 4 5 23 
HPC Cover - 12 3 7 
Stop Valve - - l ~ 
Regulator Valve - i) 2 7 
Housing 5-15 l 11 16 12 
Cover 18 8 1] 6 
Stop Valve 4 l 4 3 
Regulator Valve 12 4 20 9 
Housing 15 — 30 ~ ] 7 4 
Cover 2 - 4 l 
Stop Valve ~ - 5 6 
Regulator Valve l ~ 3 _ 
Housing 30 — — 6 4 
Cover - - - 3 
Stop Valve 2 l 3 l 
Regulator Valve 3 l l l 


An investigation of the cracks in the high-temperature region of the housing showed that the cracks are 
characteristic of heat-fatigue, are filled with oxides, and have pinched, dead ends (Fig. 1). This character 
of the cracks testifies to their thermal origin. 


We should note that the damage in each of the power stations is not identical. Thus at one power station 
there are practically no defects at all on cast parts, while conversely at another there are defects on 
practically all the cast parts. 


The most damaged area of the housing is the exterior surface of the chamber of the regulator stage. 


In recent years VTI and TsKTI have developed a method for determining the crack resistance of the 
metal in cast turbine housing parts®’*. The crack resistance characteristics were determined by this 
method for the metal of three turbine housings after 200,000 h of operation. 


2k 








oe 
a 


wee 4 


ed 
. . 4 5 
~ Bit $5 
Hen 


Puc. 8 





The mechanical properties estimated from standard samples are given below (1 is for no cracks over the 
entire period of operation, 2 is for an average number of defects — 12, and 3 is the maximum number of 
crack formations observed over a.! the 13 turbines — 25), 


l 9 8 PMTU-17-56 
a (LMZ) 
op, MPa 470.0 440.9 485.0 470.0 
o¢,MPa 314.0 300.0 268.0 250.0 
5,% 26.0 26.0 26.0 18.0 
wv, % 59.0 51.0 59.0 30.0 
ay kd /m? 700 98 147 290 
Hp 170 143 149 135 — 180 


The hot hardness H p and critical coverage o , determined by the VTI method, and the critical coverage 
0 determined by the TsKTI method are presented below. 


Hp 7«,mm To 
l 1200 0.25 0.25 
2 1000 0.40 0.20 
3 1100 0.35 0.22 


From an examination of this data it follows that the crack-resistance characteristics for 20KhML steel in 
all cases exhibit satisfactory values. 


A study of the metal after 275,000 h of operation at the recommended capacity shows that the metal is in 
satisfactory condition®. In addition to the accepted methods, the high-pressure cylinder metal was 
subjected to an evaluation of damage in terms of the criteria proposed in 67 The damage was 
determined from electron microscope investigation of foils, electron fractography of low-temperature 
breaks, and optical metallography. Viewing a large number of foil sections revealed that damage to the 
metal was low, and only a small number of individual micropores 0.1 — 0.3 zm in size were detected along 
the boundaries of the ferrite grains (Fig. 2). 


We should note that damage to the metal as estimated from the dislocation structure of ferrite was not 
taken into consideration, since the metal has a large-grain structure and there is practically no separation 
of carbides by the main part of a grain. Thus fragmentation and formation of cell-like structures were 
noted in the metal. 


Investigation reveals that the breaks are brittle and mixed, encountering both transcrystalline sections and 
sections of between-grain breaks. Only a few fine pores were found by examining the breaks on a 
scanning microscope, which established that the grain boundaries are hardly damaged at all (Fig. 3). 


Metallographic analysis on an optical microscope detected no pores. 
Thus this investigation of housing parts has shown that in terms of all the indicators — crack effects, crack 


resistance, microstructural conditions, mechanical properties, phase composition and carbide formation 
processes, in terms of damage the metal is in satisfactory condition. 
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High-Pressure Rotor. Now we will conduct a fairly detailed investigation of the struc ural condition of 
the metal in the rotor after extended operation®, which revealed an absence of substantial alteration in the 
structure and properties. We followed the condition of a large number of high-pressure rotors using non- 
destructive monitoring methods. Practically all the checked rotors were free of defects. Defects were 
detected in all on only a few rotors. Investigation of the nature of the defects showed that these cracks 
were metallurgical in origin and located near segregated material; long-term operation does not lead to 
them opening any further. 


Fig. 4 presents a typical crack 4 mm deep found in the regulator stage area, which hardly grew at all after 
210,000 h of operation. This leads us to the conclusion that the level of active stress in the rotor is 
substantially less than what has been calculated. 


Thus the stability of the properties of 34KhMI1A steel during extended operation and the absence of 
action on the rotor of a cyclic stress nature that would cause crack generation all point to the possibility 
of further reliable long-term operation of the high-pressure rotor. 


Low-Pressure Housing. When following the condition of housings after 200,000 h of operation, we noved 
considerable corrosional and erosional wear. Washings were found on the cylinder baffles, washouts of the 
guidevane roots along the inner generatrix and throat of the diaphragm, washings on the diaphragm 
joints and at the guidevane roots, and wear of the flow separator (Fig. 5). 


Low-Pressure Rotor. Washouts were found on the disks along the outside near the rivets on the steam 
intake side (Fig. 6). 


Annular cracks were located on the radius transitions from the circumference to the balance band. Radial 
cracks were found near the rivet holes (Fig. 7). Cracks were found in the keyways (Fig. 8). There were 
washouts in the hub of the disk below the diaphragm seals extending to the disk apron. 


Thus this data indicates that the operational reliability of the low-pressure cylinder falls off sharply after 
200,000 h of operation due to corrosional/erosional wear of the disk, diaphragm, baffles and flow 
separator, as well as corrosional/fatigue cracking of the disks in 13 to 19 of the stages. 


Conciusions. 


1. The high-pressure cylinder of the VK-100-90 turbine may be reliably operated for up to 500,000 h with 
vigilant monitoring after cracks begin forming in the disk slots beneath the blades. 


2. The low-pressure cylinders are not operationally reliable after 250,000 h of operation, requiring large 
outlays during capital repairs. 


3. Turbines that are to be converted to central-heating uses do not require replacement of the high- 
pressure cylinder, as has been specified in the program, but rather of the low-pressure cylinder. 
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[Text] The situation frequently arises in design of optical 
communication systems when the light guide changes from the 
high-pressure to the low-pressure zone. A sharp drop of external 
pressure, acting on the fiber, leads to intermittent variation of its 
elastic characteristics at the interface of zones, which in turn leads 
to the corresponding variation of the dielectric constant of the fiber. 
Therefore, the emission, propagated along the fiber, is reflected at the 
zone interface. It is physically obvious that energy losses of the 
information signal, caused by reflection, will be minimal if the 
pressure drop is smooth and extended on some section of the fiber. It 
should be noted that smooth variation of the stress tensor components 
contributes not only to a decrease of radiant energy losses, but also to 
a decrease of the probability that the fiber will chip on the section of 
the pressure drop, which is in itself an important problem. 


To reduce the pressure jump, a transition region of the fiber from the 
direction of the high-pressure zone must be covered by protective 
material, the elastic characteristics and configuration of which should 
be selected so that the radiant energy losses be minimal. 


The distribution of elastic stresses along the fiber on the transition 
region was studied in this paper and the conditions at which the 
pressure jump is minimal were determined. Taking the symmetry of the 
problem into account, let us write the equations of elasticity ina 
cylindrical coordinate system [1]: 
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(1) 








where v is Poisson’s coefficient of matter, u, is the radial 
displacement of the particle, and u; is the displacement of the particle 
along the axis of symmetry of the structure, coinciding with axis z. 

Let us discuss the boundary conditions to the equations of elasticity: 


a) It is obvious that radial displacements at the fiber-protective 
material interface are identical for both media: 


lis (Fy) = Ua, (ry). (2) 


Here and further subscript “1” is ascribed to the fiber, while subscript 
"2" is ascribed to the protective material, and ry is the radius of the 
fiber. 

b) Assuming that the adhesion between the fiber and protective material 


is sufficiently good, one can assume that the axial displacements at 
point r = ra for both materials are identical: 


Uy, (yh lis. UF yt. (3) 


c) The radial stresses at the interface r = re in the fiber and 
protective material are identical: 


Sir, (ly) Sup, UD. (4) 


d) The protective material-external "internal pressure” medium should 
be set equal at the interface to the external pressure [2]: 


Samy No=Pm: (5) 


where Jp, are the components of the stress tensor, ng are the components 
of the external normal to the surface of the protective layer, and pm 
are the components of external pressure. 





If the equation of the generatrix of the cone-like protective layer is 
written in the form of r = r(z), boundary condition {5) will assume the 


form 


dr 
o,,+ p=a,, ——, (o,,+p) © et. 
dz : 


ul 


(6) 


The following boundary conditions should be added to the conditions of 
(2)-(4) and (6) 


u, (r=9;=0, u,(z=Oi=0. (7) 
Let us find the solution of system (1) in the following form: 


3 


2 


DD] MN. 


i (r) + 2F (7), u,= 


a) % 


u,= 


Piri +RODir). (8) 


Substituting these expressions into equations (1), it is easy to 
determine the functions contained in formulas (8): 


: r : 
finec) In —+c), 
R 


























r 3 r r R 
glry=c) _ —~__ et) _ Jn + cl _, 
OR Tien RRIF 
| r° 3 8(1—v)?—!I | 
F (r)=— — | c%)—_ cf +3 (1--v) ch | — 
I1—2y R? | 8 ]—v 
mS cM) soe_L, ad In 1 4 cf) In $c"), (9) 
4 (I—v) R? R 
3 2 
® (r)= 7 | dye. bd — CUE SZC) | L 
1—2y 8R§ 4 |]—y 
3 5) i 
+ 3 in £ Sot sil—eve” ( In 22) 4 
l—y SR? R 4(l—v) R R 2. 
r ke 
+c) ~~, cls) 
R r r 
where cl) ((=1, 2, j=1, 2...8) are the desired constants, determined from 


the boundary conditions, and R is the cross-sectional radius of the 
protective layer at z = 0. 
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Let us further consider cases when Poisson’s coefficients of the fiber 
and protective layer are identical (v; = v2 = v). Having used the known 
formulas (see [2] for example) that link %mq to the components of the 
particle displacement vector and substituting formulas (8) and (9) into 
boundary conditions (2)-(4), (7), and (8), we find the following 
relations for the desired constants: 


cl) =o) a=nc{*) = cl) =c{ = cf) =0, cf*) ex cf?), 


2 
=e < Bid +3ve!%), chY=B (c? +3 vel), 
0 








c= + oie (c+ 3 yvci"), 
1—2y 
cl") —8 en + ve + a (4 vel +(3—4 v) 30%) | + (10) 
Bch) v(i—v) , 1—2v r, | 
—E E,{j In ——— |] — 
+E ‘129 Ti) ( "R 5) 

I—2y ry, | 

—E, ame [In —+--— ; 
‘Ty (i R 9 ’)| 


cf cf) + ne cf) — Sv ci) ( In Zs) 
r 2(1-—2v) 2 (i1—v) 2 





where the following notation is introduced 


ae 
Rt £,+E,(1—2y) (11) 





E; and Eo are Young’s modulus of the fiber and protective material, 
respectively. 


It now remains to use boundary conditions (6). It follows from general 
analysis of conditions (6) that a significant reduction of pressure on 
the fiber at z —- O may occur only with a slow decrease of function r(z). 
One can also establish that conditions (6) in the region of z = R are 
fulfilled for any z, if r(z) has the form 


2 


fan OF (2 R) (12) 


The dimensionless parameter oa in this expression determines the rate of 
decrease of the generatrix of the cone-line protective layer and, as 
follows from the above arguments, should satisfy the inequality 
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Substituting (12) into boundary conditions (6) and taking into account 
relations (10), one can calculate the constants C')’j;, and one can then 
use them to calculate the distribution along the coordinates of the 
stress tensor components in the fiber. They have the following form: 





E, ! 2a(I—2 vp 2 
Sp =Tigg =—P [ a al 
E,+(I—2yE, 1—(l—2WB I—(I—2) B R 


Qty 
E, 1—(l—2y) 6 | 


Pee [2 (1+5)-1] Bh 


(14) 


o1::;=—P 





L+2v + 
To 





9 

Terms proportional to at were not considered in derivation of formulas 
(14), since, besides small parameter a, they contain small numerical 
coefficients. 


Analysis of formulas (14) shows that a significant decrease of pressure 
at z - 0 occurs only if the parameters of the protective layer R, v and 
E» satisfy the inequality 





Thus, to reduce the energy losses of the information signal passing 
through the optical system with the section of the pressure drop, it is 
necessary that the parameters of the transition region satisfy 
conditions (13) and (15). 


BIBLIOGRAPHY 


1. Muradyan, A. G., I. S. Goldfarb and V. P. Inozemtsev, “Opticheskiye 
kabeli mnogokanalnykh liniy svyazi" [Optical Cables of Multichannel 
Communication Lines], Moscow, 1987. 


2. Landau, L. and Ye. Lifshits, “Teoriya uprugosti” [Elasticity 
Theory], Moscow, 1965. 


- END - 


33 





END OF 
FICHE 


DATE FILMED 
NOY ATU 











a sss—sssSs—‘(‘(‘(‘(Csrss—————C‘CtwrSCis 





